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A Mechanism of Spin- Triplet Superconductivity in Hubbard Model on 
Triangular Lattice: Application to UNi 2 Al 3 
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We discuss the possibility of spin-triplet superconductivity in a two-dimensional Hubbard 
model on a triangular lattice within the third-order perturbation theory. When we vary the 
symmetry in the dispersion of the bare energy band from D 2 to Dq , spin-singlet superconduc- 
tivity in the D 2 -symmetric system is suppressed and we obtain spin-triplet superconductivity 
in near the Dg-symmetric system. In this case, it is found that the vertex terms, which are 
not included in the interaction mediated by the spin fluctuation, are essential for realizing 
the spin-triplet pairing. We point out the possibility that obtained results correspond to the 
difference between the superconductivity of UM2AI3 and that of UPd2Al3. 
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tice, third-order perturbation theory, vertex correction, UIMi 2 AI 3 , UPd 2 AI 3 

Heavy-fermion superconductors have attracted much interest since spin-triplet supercon- 
ductivity has been reported in these compounds. Respective isostructual hexagonal PrN^Ahj- 
type heavy-fermion compounds, UPd2Ai3 and UNI2AI3, were found to be superconductors 
with T c ~ 2 and IK, respectively. Their superconducting states coexist with antiferromag- 
netic long-range orders with Tn = 14.3 and 4.5K, respectively. The properties of magnetism 
and superconductivity are different in each compound. UPd2Ai3 orders in an antiferromag- 
netic long-range structure with a commensurate wave vector k = (0, 0, 0.5) and has a large 
ordered magnetic moment \x = 0.85^b ° n uranium atoms. UM2AI3 orders in a spin-density 
wave (SDW) with an incommensurate wave vector k = {H ± 0.61,0,0.5) and has a tiny 
magnetic moment of 0.2/zb whose magnitude is modulated within the basal plane. Regarding 
the superconducting properties, various experiments suggest that UPd 2 Ai3 is a spin-singlet 
d-wave superconductor with a line-node gap. On the other hand, spin-triplet superconductiv- 
ity in UN12AI3 has been revealed by the NMR measurement recently. 1 We have proposed a 
theory of the superconductivity in LT^Ahj. 2 Therefore, we have to explain the mechanism of 
the spin-triplet superconductivity in UN12AI3 and the difference between the superconducting 
state of UN12AI3 and that of UPd2Al3 from the same viewpoint as that of our previous work, 
although the detailed electronic structure of UN12AI3 has not been investigated yet. Motivated 
by the point mentioned above, we discuss the possibility of spin-triplet superconductivity in 
a two-dimensional Hubbard model on a triangular lattice, as a first step. In the case of spin- 
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triplet superconductor Sr2Ru04, Nomura and Yamada have recognized that the momentum 
and frequency dependence of the effective interaction between electrons, which is not included 
in the interaction mediated by the spin fluctuation, is essential for realizing the spin-triplet 
pairing and they have explained the superconducting mechanism within the third-order per- 
turbation theory (TOPT). 3 The perturbation approach is sensitive to the dispersion of the 
bare energy band, by its nature. It implies that the lattice structures and the band filling play 
essential roles in the calculation of superconducting transition temperature T c . Therefore, it 
is important to evaluate superconducting transition temperature T c on the basis of the de- 
tailed electronic structure in each system. In this paper, we also calculate T c of spin-triplet 
superconductivity in a two-dimensional Hubbard model on a triangular lattice within the 
TOPT. In a similar model, Kuroki and Arita have proposed that spin-fluctuation-mediated 
spin-triplet superconductivity can be realized in their model with disconnected Fermi sur- 
faces, by using fluctuation-exchange approximation (FLEX). 4 However, as mentioned later, 
our proposed mechanism of spin-triplet superconductivity is completely different from their 
proposed mechanism. 

We start from the quasi-particle state according to the discussion of our previous work. 5 
We write the model Hamiltonian as follows: 

H = Hq + Hu (1) 

H = ^( e (k)-/i)4 CT a kCT , (2) 

k,cr 

Hl = 2N ^2 ^2 ^k 1 +k2,k3+k4«ki CT a k 2 a' a k3a'«k 4CT > (3) 
cr^cr' kik 2 k3k4 

where a ko .(ak<r) is the creation(annihilation) operator for the electron with momentum k and 
spin index a; e(k) and fx are the dispersion of the bare energy band on a two-dimensional 
triangular lattice and the chemical potential, respectively. The sum over k indicates taking 
summation over a primitive cell of the inverse lattice. In the above equations, we have rescaled 
the length, energy, temperature and time by a,t, j. (where a, t, fee and h are the lattice 
constant, the nearest neighbor hopping integrals, Boltzmann constant and Planck constant 
divided by 2ir, respectively) 

We calculate T c by solving Eliashberg's equation (Fig. 1). In the equation, the normal 
self-energy and the effective interaction are obtained within the third-order perturbation with 
respect to J7(Fig. 2). The diagrams enclosed by a dashed line in Fig. 2(b) are the vertex 
correction terms which are not direct contributions from spin fluctuations. The other diagrams 
are included in RPA. We call the latter 'RPA-like diagrams' in this paper. In Fig. 2(b), we 
omit writing the diagrams given by turning the vertex correction terms in Fig. 2(b) upside 
down. 

Our model parameters are the dispersion e(k) of the bare energy band on a two- 
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dimensional triangular lattice, the electron number n per one spin site and the Coulomb 
repulsion U. Regarding the dispersion of the bare energy band, we consider the following 
dispersions of the bare energy band. 
At first, we consider the following: 

\/3 1 

e D2 (k,t m ) = -4cos(— k x )cos(-k y ) - 2t m cos(k y );t m / t(= 1). (4) 

6D 2 (k, t m ) exhibits only the D2-symmetry because of t m / 1 (Fig. 3(a)). In the case of spin- 
singlet superconductor UPd2Al3, we adopted the above dispersion of the bare energy band in 
the previous work. 2 In the previous work, we have considered only the nearest neighbor hop- 
ping integrals and we have assumed that the value of hopping integral t m along the magnetic 
moment is different from the value of other hopping integrals t, because the superconductivity 
of UPd2Al3 is realized in the antiferromagnetic state. Thus we have included the effect of the 
antiferromagnetic order in the difference between t m and t. We have also determined the val- 
ues so as to reproduce the considered Fermi sheet which is obtained by the band calculation 
and is not of Dg-symmetry, reflecting the antiferromagnetic structure. We have concluded that 
the main origin of the superconductivity is the momentum dependence of the spin fluctua- 
tions which stems from the shape of our considered Fermi sheet which undergoes symmetry 
breakdown(D6 — > D2-symmetry) due to the antiferromagnetic order and then possesses nest- 
ing properties. In the present paper, we investigate how superconducting states change when 
we vary t m (symmetry in the system) from t m / 1(D2) to t m = 1(Dq). We also investigate 
the possibility of spin-triplet superconductivity near the previous model of UPd2Al3. Then 
we discuss a possible model of UM2AI3, although the detailed electronic structure of UN12AI3 
has not been investigated yet. 

Next we consider Dg-symmetric dispersions of the bare energy band (i.e, dispersions of 
the bare energy band without anisotropic nature in the hopping integrals). In this paper, we 
consider the nearest neighbor hopping model. We also consider the following dispersion of the 
bare energy band: 

e D6 (k) = e D2 (k, 1). (5) 

To satisfy Luttinger's theorem, that is, the conservation law of particle number, we adjust 
the chemical potential fi by using the secant method. To solve Eliahberg's equation by using 
the power method algorithm, we have to calculate the summation over the momentum and the 
frequency space. Since all summations are in the convolution forms, we can carry them out by 
using the algorithm of the fast Fourier transformation. For the frequency, irrespective of the 
temperature, we have 1024 Matsubara frequencies. Therefore, we calculate throughout in the 
temperature region T > T\ im , where T\ im is the lower limit temperature for reliable numerical 
calculation, which is estimated to be about 3.0 x 10~ 3 (> Ae/(27r x 1024) ~ 1.4 x 10~ 3 ), where 
Ae is the bandwidth; we divide a primitive cell into 128x128 meshes. 



3/?? 



J. Phys. Soc. Jpn. 



Full Paper 



Table I. Irreducible representations of D 2 



Irreducible representation Parity Basis functions (maximum wavelength in k-space) 

A even 1 

-a (!) (!) 

Bi even c\ — c 2 

b 2 o^ S S 1) + 4 1) ,4 1) 

B 3 odd s ( l ] - 



To investigate how superconducting states change when we vary t m (symmetry in the 
system) from t m / 1(D2) to t m = l(De), we calculate the eigenvalues of Eliahberg's equa- 
tion for various i m (0.75 ~ 0.99) and n, starting from the model parameters of our previous 
work (eo 2 (k, t m ), t m = 0.75, n = 0.572). In this calculation, we select the Coulomb repulsion 
U = 7.5. In the case of the D2-symmetric system, we can classify the eigenvalues of Eliahberg's 
equation according to the irreducible representations of D2 . D2 has four irreducible represen- 
tations (see Table I). 

The short notations and in Tables I and II mean that = sin(R^ • k) and 
Cj^ = cos(R^ • k) , where is defined in Fig. 3(b). 

The eigenvalues which belong to A or Bi correspond to spin-singlet states. The eigenval- 
ues which belong to B2 or B3 correspond to spin-triplet states. Based on the values of our 
parameters, in the case of spin-singlet states, the maximum eigenvalue belongs to Bi and in 
the case of spin-triplet states, the maximum eigenvalue belongs to B2. 

The results are shown in Fig. 4. From this figure, the following facts are evident. 

In the case of n = 0.572 and t m = 0.75, the spin-singlet state is most stable. We have 
pointed out that this spin-singlet state corresponds to that of UPd2Ai3 in our previous work. 
When we vary t m (symmetry in the system) from t m = 0.75(D2) to t m = l(De), at the same 
time, the anisotropic nature in spin fluctuation is suppressed. In this case, the spin-singlet 
state is suppressed because the main origin of the d-wave superconductivity is the momentum 
and frequency dependence of spin fluctuations as we have pointed out in our previous work, 
and we can see that spin-triplet states have the tendency to emerge. 

Next, we consider Dg-symmetric dispersions £D 6 (k) of the bare energy band and calculate 
T c for some model parameters U and n. In the case of the Dg-symmetric system, we can 
classify the eigenvalues of Eliahberg's equation according to the irreducible representations of 
D6- Dg has six irreducible representations (see Table II). 

The eigenvalues which belong to Ai, A2 or E2 correspond to spin-singlet states. The 
eigenvalues which belong to Bi, B2 or Ei correspond to spin-triplet states. Based on the 
values of our parameters, the maximum eigenvalue belongs to Ei which corresponds to a 
spin-triplet state. 
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Tabic II. Irreducible representations of D 6 



Irreducible representation 


Parity 


Basis functions (maximum wavelength in k-space) 


Ai 


even 


1 


A 2 


even 


J4) .(4) (4) (4) (4) (4) 
L l L 2 3 4 c 5 L 6 


E2 


even 


{c<"-4 1) ,4" + c<"-2c<"} 


Bi 


odd 


4" + 4" +4" 


B 2 


odd 


4 2) +4 2) +4 2) 


Ei 


odd 


{.<•'- 4", 4" +4" -24"} 



The dependences of T c on U and n are shown Fig. ??. On the basis of these results, we can 
point out the following findings. For large U, high T c values are obtained for all parameters. 

When we increase the electron density n for the fixed valued Coulomb repulsion U and 
the Fermi level becomes close to the van Hove singularity (whose results are not presented in 
this paper), then high T c values are obtained. 

To examine how the vertex corrections influence T c , we also tried to calculate T c by 
including only the RPA-like diagrams but we could not find any finite value of T c within the 
precision of our numerical calculations. Therefore, we can calculate only the eigenvalue by 
including only the RPA-like diagrams at T c obtained by using the TOPT. 

From these results, it is found that the vertex terms, which are not included in the inter- 
action mediated by the spin fluctuation, are essential for realizing the spin-triplet pairing. 

In this paper, we have discussed the possibility of spin-triplet superconductivity in a two- 
dimensional Hubbard model on a triangular lattice within the TOPT and we have obtained 
spin-triplet superconducting states in near a Dg-symmetric system. In this case, the vertex 
correction terms play an essential role in obtaining the spin-triplet superconducting states. 
This is the same result that Nomura and Yamada have obtained in the investigation of the 
mechanism of spin-triplet superconductor Sr2Ru04 within the TOPT, as we have mentioned 
above. Jujo et al. have investigated the mechanism of the spin-singlet superconductivities 
in organic superconductors K-type (BEDT-TTF) 2 X within the TOPT and FLEX. 6 They 
have compared the results between the TOPT and FLEX, and have pointed out that the 
vertex correction terms (which are also not included in the FLEX) have a crucial effect on 
the calculation of T c for strongly frustrated systems. Therefore, our presented results are 
consistent with the above two results. Recently, Nisikawa et al. have discussed the spin-singlet 
superconductivity in Celr-cCoi-^Ins within the TOPT and have pointed out that the vertex 
correction terms are important to explain the x-dependence of T c . 5 The findings described 
above suggest that the calculations of T c , which include only the spin fluctuation terms, are 
questionable and should be carefully performed with vertex corrections. 
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UNi2Al3 and UPCI2AI3 undergo superconducting transition below the Neel temperature. 
Therefore, the symmetry in the itinerant electron system under the antiferromagnetic struc- 
ture is one of the important matters to consider in investigating the mechanism of super- 
conductivity. In the case of UPd2Al3, the symmetry in the itinerant electron system is not 
of hexagonal symmetry, reflecting the effect of the antiferromagnetic structure with a large 
ordered magnetic moment jjl = 0.85^b on uranium atoms, and we have treated this by consid- 
ering the anisotropic hopping integral t m 7^ 1. In this context, the symmetry in the itinerant 
electron system of UM2AI3 may be more isotropic than the symmetry in that of UPd2Al3 
because of a reflecting incommensurate SDW order with a tiny moment of 0.2^b 5 although the 
detailed electronic structure of UM2AI3 has not been investigated yet. Therefore, we assume 
that t m ~ 1 in the case of UN12AI3. Based on the hypothesis mentioned above, our results 
seem to explain not only the mechanism of spin-triplet superconductivity in UN12AI3 but also 
the difference between the superconductivity of UM2AI3 and that of UPd2Ai3 because the 
spin-singlet superconductivity in the D2-symmetric system(t m 7^ 1) is suppressed toward De- 
symmetry and the spin-triplet superconductivity emerges in near the Dg-symmetric system 



In conclusion, we discussed the possibility of spin-triplet superconductivity in a two- 
dimensional Hubbard model on a triangular lattice within the TOPT. We obtained spin-triplet 
superconducting states in near the Dg-symmetric system. We pointed out the possibility that 
our obtained results correspond to the difference between the superconductivity of UN12AI3 
and that of UPd 2 Al 3 . 

Numerical computation in this work was carried out at the Yukawa Institute Computer 
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Facility. 
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Fig. 1. Eliashberg's equation. The thick line represents Green's function with self-energy correction. 
The shaded rectangle represents the effective interaction. 




Fig. 2. (a) Feynman diagrams of the normal self-energy up to the third order, (b) Feynman diagrams 
of the effective interaction up to the third order. Solid and dashed lines correspond to the bare 
Green's function and the interaction, respectively. 
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Fig. 4. (a) The calculated T c is shown as a function of t m for various electron numbers n per one 
spin site, (b) The calculated eigenvalue is shown as a function of t m at Ti; m = 0.003 for n = 0.572 
and 0.600 
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Fig. 5. The calculated T c is shown as a function of U for various electron numbers n per one spin 
site." In this figure, (v=***:RPA-like)" means that eigenvalue v calculated by including only the 
RPA-likc diagrams is v =*** at T c obtained by using the TOPT. 
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